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Abstract 

We calculate the leading divergences at NNLO for the octet part of the nonleptonic 
weak sector of chiral perturbation theory, using renormalization group methods. The 
role of counterterms which vanish at the equation of motion and their use to simplify 
the calculation is shown explicitly. The obtained counterterm Lagrangian can be 
employed to calculate the chiral double log contributions of quantities in this sector, 
most notably the K — > tttt amplitude. The double log contribution of the latter is 
discussed in a separate paper. 
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1 Introduction 



We determine the leading divergences at NNLO for the nonleptonic weak chiral Lagrangian 
which transforms as an octet under the chiral group, extending the NLO calculation of the 
latter |2] . The obtained counterterm Lagrangian can be used to calculate the leading 
logarithmic contributions, double logs in short, of observables in this sector. These con- 
tributions are in particular interesting since all the low energy constants of higher order 
are unknown in this sector. The calculation of logarithmic contributions provides thus the 
only way to get a first analytical estimate of the size of the higher order corrections one 
has to expect. Analogue calculations in the strong sector of chiral perturbation theory 
(CHPT ) have already been worked out j^Ej. The results of this paper are used to calcu- 
late the double logs to the K — > irir amplitude, presented in a separate paper [5]. For the 
latter amplitude, methods have also been worked out to extract the needed next to leading 
order (NLO) low energy constants (LEC's) by lattice simulations; however, the proposed 
approach is rather ambitious [SJ Ej • Interesting further applications of the results obtained 
here are for instance the calculation of the double logs of the K — > nun or the K —* 7777 
amplitude. 

The chiral logs are introduced during the process of renormalization jHj. These logarithms, 
which correspond to the infrared singularities when the masses of the theory approach zero, 
can produce sizable contributions to observables. Using dimensional regularization, it is 
straightforward to understand how the leading logs are related to the leading counterterms: 
while renormalizing the theory, one has to introduce an energy scale \x to ensure the correct 
dimensions of observables calculated in d— dimensional space-time. In particular, for the 
divergences generated by loop calculations, this means that they can only show up in the 
following structures: 

To illustrate the order of magnitudes of these chiral corrections, Tabled displays the various 
contributions up to NNLO for the nn scattering lengths in two flavor CHPT , showing that 
the double log contribution in this case amounts to almost the full NNLO corrections, 
corresponding to close to 10% of the total result 



LP NLO NNLO DOUBLE LOGS 

a° 0.156 0.044 0.017 0.013 
cft-al 0.201 0.042 0.016 0.012 



Table 1: The leading order nn scattering lengths in two flavor CHPT with the chiral cor- 
rections up to NNLO, for some standard values of the NLO LEC's and at a renormalization 
scale fi = lGeV. The double logs are included in the NNLO correction. 



2 



For the case of three flavor CHPT , we show the chiral corrections up to NNLO to the 
pion and Kaon decay constants and the vector form factor of fOJ El E| i n Table El 
One notices that the relative size of the double logs is less pronounced than for two flavor 
CHPT . Typically, the double logs amount to 20 — 35% of the total NNLO contributions, 
corresponding to about 10% of the total corrections to the leading order result. Although 
the numerical values of these corrections are not too large, one should keep in mind that 
for applications like chiral extrapolations, the relative size of the double log contributions 
to the NLO corrections is of importance, which is in the range of 10% and therefore sizable. 
Let us in this context again emphasize that while it is certainly appropriate to perform the 
full two loop calculations in the strong sector, due to the lack of knowledge of the weak 
LEC's the calculation of the double log contributions is presumably the best one can ever 
achieve to get an estimate of the size of the NNLO corrections in the weak sector. 

LP NLO NNLO DOUBLE LOGS 

F n /F 1 0.068 -0.172 -0.050 

F K /F n 1 0.216 0.035 0.06 

f+(0)[K a ] 1 -0.023 0.015 0.004 

Table 2: The chiral corrections up to NNLO for the pion and Kaon decay constants and 
the vector form factor of Ki 3 . The values can however vary considerably, depending on the 
LEC's one employs. The numbers above are calculated with some standard values of the 
NLO LEC's and all the renormalized NNLO LEC's set to zero at y, = 770MeV. 

The outline of the paper is as follows: 

In section |2] we introduce the notation used and give the needed CHPT Lagrangians. Sec- 
tion El provides a very brief overview of the general framework in which this calculation was 
performed; the generating functional is introduced, and it is shown how one can use the 
background field method to calculate the counterterms needed to renormalize the latter. 
In section 0] we discuss the role of operators which vanish at the solution of the equation 
of motion (EOM terms) to simplify the calculation of counterterm Lagrangians: One can 
choose the coefficients of these EOM terms in a way that the sum of all one particle re- 
ducible (1PR) topologies at a given ft-order won't generate any divergences. We pin down 
these coefficients and henceforth only need to take into account one particle irreducible 
(1PI) topologies. 

In section we sketch very briefly the renormalization group techniques which are em- 
ployed in the present calculation. The basic result of this section is that one can obtain the 
border 2 highest pole counterterm (NNLO) by performing a one loop calculation which 
uses the border 1 (NLO) counterterm as input. In the last part, section El we illustrate 
how the concrete calculation works with two simple examples. 
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2 CHPT Lagrangian 

The lowest order chiral Lagrangian which allows for AS — 1 strangeness changing interac- 
tions is given by (Throughout this section we will work in euclidean space-time): 



L® = + £& =1 



(2.1) 



which encodes the dynamics of the pseudo-Goldstone bosons in the presence of external 
source fields s,p, v^, a M . The first term corresponds to the strong interaction Lagrangian: 



4 0) = ^(<«,«M> "CO) > 

where (■) stands for the flavor trace. Further we used: 

= i(u t (9 M -ir M )«-«(9 M -iZ M )« t ) , 



(2.2) 



u, 



[2.3) 



A list of additional building blocks used for the Lagrangians of higher order can be found 
in appendix ID1 

The u matrix encodes the octet of the light pseudo-scalar bosons in the exponential 
p ar amet r iz at ion : 



/ 7T° I ys 

1 Vz^~ Ve 



u = exp( 



V2F' 



V 



7T 

K~ 



7T+ K + \ 

A V6 ) 



(2.4) 



The definitions in Eq. ()2.3|) embody also r M and / M , the external vector source fields, 
whereas the scalar counterparts are encoded in \ iy^ being the quark mass matrix): 

X = 2B (^ + s{x)+ip(x)) . 

Bq is related to the vacuum expectation value of the scalar quark density: 

(0\qq\0) = -F^B (l + O(^)) . 

Fq ~ 92.4MeV corresponds to the chiral limit value of the pion decay constant. 
The Lagrangian which triggers flavor changing processes is given by: 



L 



(0) 

AS=1 



CF^(g 8 (A 32 u^) -g s (A 32X+ ) + g 27 t^ kl (A ll u fl )(A kl u ll ) + h.c. , (2.5) 



with t 11 ' 23 = t 13 ^ 21 = t 21 ; 13 = t 23 ^ 1 = 1/3, t 



22:23 



t 



23;22 



t 



23;33 



t 



33;23 



■1/6, and all 
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other t's vanishing. 
The constant C: 




(2.6) 



Aij := uXijVj ; (Xij)ab = SaiSbj 



The first two operators in Eq. (|2.5|L proportional to g$ and g' 8 transform as an octet, (8, 1), 
under the chiral group SU(3)^ <S> SU(3)^, whereas the last, proportional to g 2 7, transforms 
like a 27-plet, (27,1). We neglected a third contribution transforming like (8,8), which 
takes into account virtual photons. 

Since the octet part of Eq. ()2.5|) is believed to be the main source of the AI = 1/2 rule, we 
will only use the latter in our calculation. Furthermore we discard the part proportional 
to g 8 , since it doesn't contribute to on-shell processes [T2*l IT3*1 ITU Q]. Of the remaining 
operator we only use its CP even part, which we will henceforth denote by 



In addition to the lowest order Lagrangians £j s and Lw discussed above, we will also use 
the NLO Lagrangian Lw\ introduced in Eq. (|3.6)l . 

3 The h expansion of the generating functional 

In this section we set up some notation and discuss the NLO and NNLO expressions for 
the generating functional of the nonleptonic weak chiral Lagrangian. 

3.1 Notation 

The generating functional is defined as the vacuum to vacuum transition amplitude in the 
presence of sources, collectively denoted by j (Throughout this section we will be working 
in euclidean space-time.): 



:= CF£gs(Au^) 



A := u\qu^ . 



(2.7) 




(3.1) 




(3.2) 



oo 




(3.3) 



n=0 



n 
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and similarly for the weak part. 

The tree level generating functional Z^ [j] corresponds to the action with the lowest order 
chiral Lagrangian, evaluated at the EOM: 

Z(%) = S<®[<p,j] . (3.4) 

One can calculate the separate fo-order contributions to the generating functional by use 
of the background field method, which is briefly outlined in appendix [XJ The strong and 
weak nonleptonic chiral actions are expanded in quantum fluctuation fields £ around a field 
(f, which in the following is assumed to be a solution of the classical equation of motion, 
i.e. Sq = 0. To keep the notation simple, we will henceforth suppress the arguments of 
z [j], S[tpJ] an d 

(2) 

This expansion provides the vertices which will be needed for the calculation of Zw , Eq. 
()3.9j) . The Latin indices i, j, .. correspond to an SU (N) index as well as a space-time degree 
of freedom. The field <p is treated as a background field, and the £-field is employed as new 
integration variable in the path integral, Eq. (J3.1|) . The perturbative evaluation of this 
"new" generating functional results in vacuum diagrams with respect to the £ fields. 

3.2 NLO: 

The counterterms needed to renormalize Z^: 

Z« = ififfGy + S? , (3.5) 

were first calculated by Kambor, Missimer and Wyler p]. Throughout this paper, we will 
however use the basis of operators given by Ecker, Kambor and Wyler (EKW) 0, who 
used the EOM to reduce the former. This Lagrangian assumes the form: 

37 

^ = CgsF^c^wU , (3.6) 

i=l 

where the bare LEC's cp are split into a renormalized and counterterm part: 

cP = ^cr(c^(^e) + a[ 1) l A) . (3.7) 

is the renormalization scale and the constant c parametrizes the regularization prescrip- 
tion ( ln(c) = exp(— (ln(47r) + T (1) + l)/2 for MS ).) In addition we use the notation: 

e:=A-d; N := (4vr)- 2 ; A := — . (3.8) 

e 
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-1/4 + 7V/4 
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-N/4 


-3/4 
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10 


(A% 2 + ) 


-3/N + N/A 
-1/2 -2/iV 2 


-1/4 
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(Ax + >(x + > 
(Ax 2 ) 


-13/18 




12 








* 


13 


(Ax_XX-> 








* 


23 


(A M { X _,^}) 








* 


36 


(A[x + ,X_]) 


-l/JV + JV/4 


5/12 





Table 3: List of operators needed for the K — > nix amplitude as well as those which can be 
shifted by EOM terms, marked by an asterisk. The a^J are given for Minkowski space-time. 
N corresponds to the number of flavors. 



In Table EJ we list the operators of the EKW basis needed for the K — > -kit amplitude 
at NLO, as well as those which can be shifted by terms which vanish at the equation of 
motion (marked by an asterisk): The operators given above differ slightly from the original 
EKW basis, which used W^q = (A([x + ,xJ + x\ ~ X 2 _))- We use the above definition of 
W^q since it simplifies the discussion about the operators which vanish at the solution of 
the equation of motion in section 0] somewhat. 



3.3 NNLO: zS ] 

At border 2, we have the following diagrams: 

Z^ = —-vo^ k Gi r Gj S GktSQ St + -Gij'(v^ kl Gki + tt^/GV,- — -xv^GijGkiSi 

o 2 2 

1 1 1 

— -LrijW Q (jrfc r S ^st T ~^ m < - Jr ifc < ~ r j« s Or mn Sg U rs — — XV^irSQ Ur s t 



-0(G%) , (3.9) 



where Gij is the propagator corresponding to the £ field, whose ultraviolet divergent part is 
provided in appendix El The subscript of the vertices denotes their border . The diagrams 
corresponding to Eq. (J3.9)) are drawn in Fig. ^ 

The NNLO Lagrangian Lw\ represented by diagram k in Fig. ^ has to cancel the diver- 
gences which are generated from the loop part of Z^- It takes the form: 

*S> = ^ 8 £cfwf , (3.10) 
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j k 

(2) 

Figure 1: Diagrams contributing to the generating functional Zw ■ They split into the 
class of 1PI diagrams (a,b,d,g), the 1PR diagrams (c,e,f,h,i,j), and the fr-order 2 action S™, 
diagram k. 

with the bare LEC's: 

c? = (^(cf r (^e)+a^(^^e))A + a^A^ . (3.11) 

(2) 

3.3.1 The connection between A 2 and the double chiral logs 

The highest pole of £/ 2 \ := Ea^Wf 3 , can be used to calculate the double chiral 
logs which are generated from genuine two loop diagrams of a process under consideration. 
Let us outline how this works: The sum of all diagrams in Fig. ^with the exception of the 
counterterm diagram k will result in an expression which is proportional to the square of 
Q, Eq. ()1.1|) plus other contributions which are not related to double logs (abbreviated by 
the dots): 

c*£ajvf ) a 2 + ... = c^ . 

i i /i /i 
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The A 2 divergences above have to be canceled by A 2 , which translates into the following 
identity: 

c w - s E*'" (f M^)) 2 = -c 98 /.-*£4M 2 '(4iog(^)) 2 , 

= -c^-^fd^ff )) 2 • 

In addition to the "genuine" double logs above, there are also contributions from one 
particle reducible topologies, the LSZ reduction, plus shifts of bare parameters to their 
renormalized values in lower order contributions of the process under consideration. 



4 Equation of motion terms 

In [13] it was shown that one has to allow for operators which vanish at the solution of the 
equation of motion (EOM terms in short) to define a basis in which one-particle-reducible 
(1PR) graphs contributing to the generating functional do not generate divergences. The 
equation of motion in euclidean space-time reads: 

X := V M u M +^L = , (4.1) 

with: 

:= dpUp + [r M , Ufj\ , 

T M := - (V (<^ - ir>) + uidp - w f ) , 
X- ■= X- - (X-)/N . 

Before we can start with the actual calculation discussed in section |3] and |Hl we need to 
define this proper basis, relevant for our computation: 
Defining the generating functional of proper vertices: 

r<f)[j] = itV^W : = zi n)1P V] ;a = s,w, 

evaluated at the solution of the EOM, and its functional derivative 

rW*[J] = (rWV^])^ ; a = s,w , 
we can decompose zffi into a 1PI and 1PR part: 

4 2) = 4 2)1PI + 4 2)1PR = fL 2) - f W'Gyff*' . (4.2) 

As Eq. (|4.2j) illustrates, the 1PR portion of Z will not contribute to divergences if in 
addition to T itself all its functional derivatives are finite. The latter can be achieved by 
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appropriate additions of EOM terms to the Lagrangian [T3] (We will however conclude this 
section with a stronger proposition on this point) . 
In the strong sector, we have at fi-order 1: 

£« := £ GL + x [ ?(x_X) + xP(XX) , (4.3) 

with £ GL the usual Gasser-Leutwyler Lagrangian ^1 HH and two additional EOM terms 
with coefficients x[ , x£ ■ 

It was shown in [3| that the Is* is finite if one discards the EOM terms altogether (i.e. 
Xi = X2 — 0), which is sufficient for our purposes. This result is to be expected since 
the only building block in £ GL which corresponds to EOM terms, has a vanishing 

divergent counterterm. 

In the weak sector, things get a little bit more involved. In addition to the EKW Lagrangian 
given in Eq. (jH.fi)) . we have six EOM terms: 

W = Cg*F* ( £ c?> Wf j + E W) , (4-4) 

i=l i=l 

listed in Table EJ 



i jjg fc (iy fc (1) ) 

1 i(A[w 2 ,X]) 9 

2 i(A{ X _,X}) 12 

3 (All) 12 

4 i(AX)( X _) 13 

5 (A M K,X}) 23 

6 i(A[ X+ ,X}) 36 



Table 4: All EOM terms for £jff • In the last column we provide the operators of the EKW 
basis, given in section |3J Table HU which corresponds to the respective EOM term. 

In order to pin down the coefficients , we calculate the functional derivative T$ % = Z$ % 
explicitly: : 

= G jk - UtG^G^ + ft*. (4.5) 

The diagrammatic representation of Eq. ()4.5j) is shown in Fig. |21 To compute the ultra- 
violet divergent part of these diagrams, we have to expand the lowest order Lagrangians, 
Eqs. (|2.2j) and ()2.7j) . together, with the counterterm Lagrangian A\ = ^ZiWi, Eq. I3.fil 
In appendix |E] we give the expansion of the building blocks defined in ID. II in terms of the 



10 




Figure 2: Graphical representation of Eq. 14.51 
quantum fluctuations £ and the background fields. 

The Wick contractions are then performed using the heat kernel representation of prop- 
agators, briefly explained in appendix El For the tadpole diagram tt>^ k Gjk, one uses the 
identities given in Eq. (|B.8[) . whereas for the diagram Sq GjiGkmR) 1 ™, one employs the 
identities for products of propagators, provided in Eqs. (|B.5[) . ([B.6[) ) and (|B.9|) . A more 
explicit discussion of how the computation works is given in section EJ 
The calculation, whose result is too lengthy to be displayed here, yields for the coefficients 
e- 1 ^ ( euclidean space-time, N is the number of flavors): 

e? = j I e« = e<" = «£> = 4" = ; e<" = + | . (4.6) 

This shift of the original EKW Lagrangian corresponds to a replacement of the following 
building blocks: 

W 9 = (A[ X _,u 2 ]) — > ^-2i:=(A[V^,n 2 ])-2i , 
^ 36 = (A[X + ,X_]) - W^-2i:=(A[ X+ ,V^]>-2i . (4.7) 

It is striking that Wg and W^q are the only operators which can be shifted by EOM terms 
and have a non- vanishing counterterm coefficient ap (see Table OJ). This observation leads 
to the conjecture that V^it^ instead of x_ should be used in loop calculations. A more 
general discussion of this point will be given in a separate paper 



5 The calculation of A 2 

For the computation of A\ , we use renormalization group techniques: If we write L w , 
given in Eq. (|3.1(J|) . split into the renormalized and counterterm part as follows: 

^ 2) = ^ 8 (^ 2)r +^ 2) A + 4 2) A 2 ) , 
the RGE imply the identity: (d| n) = <9/<9cJ n)r ) : 

A? = ~a«£>«4 2) , (5-1) 
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with a\ being the counterterm of £j w defined in Eq. f!3.7jl . For the case of the weak 
nonleptonic sector, Fig. ^ shows all diagrams which can contribute to divergences at h- 
order 2. If we act with the operator on the RHS of Eq. (j5.1|) . a^d^\ on all these diagrams 
and neglect the 1PR topologies (see section HJ), there is only diagram d which can give a 
no n- vanishing contribution. The general diagrammatic representation of this statement, 
Eq. (J5.1j) . is shown in Fig. El A thorough discussion and derivation of this RGE approach 
can be found in [T5|. The specific result, Eq. (|5.1j) . was already used in |19j . 
Due to the RGE, we can therefore obtain the leading poles at /i-order 2 by computing only 
the one loop diagram d in Fig. Q]and weighting it with the factor 1/2, instead of having 
to compute the genuine 1PI two loop diagrams a and b. 




Figure 3: Calculation of A 2 with the RGE. 

In appendix [U] we provide all the operators of A^ which can contribute to amplitudes 
involving at most four pseudo-Goldstone particles and no vector- sources, denoted by A 2 . 
This result has been used for the calculation of the double log contribution to the K —>■ 
nn amplitude, which is presented in a separate paper The full expression is approx- 
imately four times the size of the truncated one shown in the appendix, and thus too 
lengthy to be displayed in this paper. It can however be obtained from the author. An 
overview about the workings of the actual computation can be found in section El 



6 Outline of the computation 

In this section we provide a sketch of how the whole calculational machinery used in this 
paper works. 

6.1 Tadpole graphs 

(2) 

We will illustrate the computation of a tadpole graph by computing a part of A 2 with 
the help of Eq. ()5.1|) . Our starting point is the counterterm of the NLO Lagrangian L^w 1 
PUI2]. For completeness we provide the expanded building blocks needed in appendix lEl 
We decided, however, not to reproduce the whole expanded expression of A± here, since 
it is rather long. 

We will rather restrict ourselves to the first building block which occurs in A[ 1 , W\ = 
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(Au 2 u 2 ): We expand W\ in the quantum fluctuation fields £: 
(Au 2 u 2 ) = (Au 2 u 2 ) 

+ 1 -(A[u 2 u 2 , x^e - (A{u 2 , {u„ A 4 }})e; t 

-i(A({u 2 , {u„ [[« MI Ai], Aj]}} + [[uV, Ai], A,])}^ 

+(A({« 2 , AjAj} + {u MJ A,}K, A,}))^e + 0(e 3 ) • (6.1) 

Due to delta functions generated by functional differentiation, no space-time integration 
survives and all building blocks and ^-fields in Eq. (J6.1|) are evaluated at the same space- 
time point x. Indices in SU(N) space are denoted by The ^-fields of the bilinear terms 
proportional to , £ l H and are then contracted with the heat kernel representations 
of Ga(x,x), d^GA(x,y)\ v=x and d^d^G^x, y)\ v=x respectively, listed in Eq. (|B.8|) . For the 
first contributing term in Eq. (|6.1j) we get the following contribution to the action S™' 



d d x(A({u 2 , {u„ [[«„, Ai], Xj]}} + [[u 2 u 2 , A*], Xj]))G A [X,X)ij 
(c/i)~ £ A 



4 

(c/i)~ £ A 



rf d x(A({M 2 , A,], Aj]}} + [[u 2 u 2 , Ai], Aj]))(af )ij + finite terms 

^x(A({w 2 , {^, [[t2 Mj A,], A,]}} + [[u 2 u 2 , XIX,])) ■ 



4 

-^([« M ,Ai][« M ,A i ] + {Ai,A i } X+ )+f.t. . (6.2) 

For the subsequent contraction of the SU(N) indices i,j one uses the completeness rela- 
tions: 

N 2 -l 

^2 (XiAXiB) = -—(AB) + 2(A) (B) , 

i=l 
N 2 ~l 

(\A)(\B) = 2(AB) - -(A)(B) . (6.3) 

This last step of the computation is obviously straightforward, and we forbear to display 
the final result, since it is again rather lengthy. 



6.2 Beyond the tadpole 

The computation of the functional derivative of Tw\ Eq. (|4.5|) . involves a diagram with 
two propagators and is therefore a little bit more involved: in addition to the SU(N) 
contractions one has to deal with the space-time dependent part of the product of the two 
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propagators. Let us again restrict ourselves to the simple case of a part of the computation 
where the vertices do not carry any derivatives acting on the ^-fields. The structure of 
such a piece is then: 

Q°= [ d d xd d yv a J k (x)G jl (x,y)G km (x,y)v l ^(y) ■ (6-4) 



Here i,j,k,l,m are again pure SU(N) indices, and a corresponds to the space-time point 
Xi as well as to the SU(N) index %. The space-time dependent part can be evaluated with 
the help of Eqs. ()B.6j) and ()B.9j) (We suppress the vertices v a J k {x) and v l ™(y) in this step): 

d d xd d yG jl (x,y)G km (x,y) = ( ^^Z W^ f ' 

d d xd d yaQ(x,y)jiaQ(x,y) km \x - y\~ d+£ + f.t. , 

IT-?/ 2 ) 1 ' T(d/2) 



d d xd d ya£(x,y)jia£(x,y) krn n e 5 d {x - y) + f.t. , 
= 2(c/i)- £ A j d d x5 3l 5 km + Lt. . (6.5) 
Reinserting the vertices again, we finally get (v aj (x) = 5 d (x — Xi)v i: ' k (x)): 



Q a = Q\xi) = 2(c/i)- £ A J d d xv a s jk (x)vi k {x)+Lt. , 

= 2{c fJl )- £ Av i J k {x i )vi k {x i )+lt. . (6.6) 

If ^-fields with derivatives are contracted, one will in addition generate derivatives acting 
on the delta- function in Eq. (J6.5|) . These derivatives can be shifted to the vertices and 
Seeley-DeWitt coefficients by partial integration. As a last step, we will again have to 
contract the SU(N) indices of the vertices with the one's of Seeley-DeWitt coefficients of 
the expansion of the propagator, Eq. ()B.5|) . using once more the completeness relations in 
Eq. (Q . 



6.3 Verification of the calculation 

This computation was exclusively performed with FORM 3.1 [20], a symbolic manipulation 
program. Since the whole calculation is thus fully automatized, one can conveniently adapt 
the code to problems whose solutions are known: In order to check the written code, we used 
it to recalculate two known counterterm Lagrangians. We replaced the original Lagrangian 
Ai (see Fig. |SJ) with the respective Lagrangians required for their calculation as starting 
point, but left the rest of the code unchanged: 
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1. We recalculated the counterterms for £«/ using the method outlined in this paper, 
i.e. by the computation of the ring diagram ~tt> u Gy instead of using the logarithm: 
|Tr(ln(A s +A„,)) and projecting out the part linear in Gf, as employed in the original 
calculation £Q. We found total agreement. 

2. We recalculated the leading poles at h-order 2 in the strong sector and compared our 
result with the one obtained by Bijnens, Colangelo and Ecker jHlll]. The outcome of 
our computation matched completely with their result. 

Since this is the very first NNLO calculation of CHPT in the weak sector, it was not possible 
to compare it directly with genuine two loop calculations. However, in our opinion, the two 
checks listed above, in particular the second one, are highly nontrivial, and yield sufficient 
evidence that the FORM code written for the computation is correct. 

7 Conclusions 

We have determined the leading divergences for the weak nonleptonic chiral Lagrangian at 
NNLO for the part which transforms like an octet under the chiral group, extending the 
analogue computation of the NLO counterterms |2j. The obtained result can be used 
to calculate double log contributions of observables, which at two loop order are the only 
quantities that do not depend on any LEC's 1 . Unlike in the strong sector, in the weak 
sector it is extremely difficult to determine these LEC's, and presumably one will not be 
able to pin them down in the near future, if ever. Thus, the double logs provide a first es- 
timate about the NNLO corrections one has to expect, without the need of these unknown 
LEC's as input. Typically, the double log contributions in three flavor CHPT amount to 
around 10% of the corrections to the lowest order result. 

Corrections to lowest order CHPT quantities are used for chiral extrapolations of lattice 
data. In these days, lattice simulations have entered a stage where one uses fully un- 
quenched quarks, and aim to predict observables with an accuracy in the range of some 
percent. In view of such high precision, it certainly makes sense to include NNLO correc- 
tions in these extrapolations. The main objective of the present calculation is the use of 
the counterterm Lagrangian for the computation of the corresponding double log contri- 
butions to the K — ► 7T7T amplitude in the A J =1/2 channel, which is presented in another 
paper [Hj. Other interesting applications are analogue calculations for the K — > mrir and 
K — > 7T77 amplitudes. 
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1 besides subtracted loop integrals which would require a fully- fledged two loop calculation. 
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A Background field method 



In this appendix we provide a very brief outline of the background field formalism [2*T] . 
which was used in the present calculation. The field if in U = exp(i\/20/-F) is split into a 
background part 0, which is normally taken to be at the solution of the equation of motion, 
and a quantum fluctuation field £: 

= + £/V2 . (A.l) 

will be used as an external field, i.e. it will not propagate; therefore, only £ can generate 
H corrections. With this substitution, the action assumes the form: 

Sty] _> S ty] + ^SWtej + ^SWSitjCk + ls[^ kl UMi+- ■ (A.2) 

In the equation above we have assumed that is a solution of the equation of motion, so 
that the linear term in the expansion vanishes. The RHS of (jA.2j) can now be viewed as 
the new Lagrangian where the new integration variables of the path integral are £ instead 
of tp. The inverse of the bilinear operator in £, S[<f)]^(x, y) = 5 d (x — y)Ay corresponds as 
usual to the propagator. Please note that we will only use the bilinear part of the strong 
chiral Lagrangian to define the propagator; the expanded form of the weak nonleptonic 
chiral Lagrangian will only be needed for the definition of vertices, each insertion thereof 
corresponding to a factor G ¥ . A is brought into the canonical form of an elliptical operator: 

Ay = (-d£ + <7(x)) - , (A.3) 

which have the following explicit form for the CHPT Lagrangian, Eq. ([2.2)1 : 

7/x ij = -gW^i'^iD ' 



1 



a 



ij 



-([up, \i\[u»,\j] + {A i ,A i }x+) • (A.5) 



and will be used to define the propagator in appendix [BJ The field strength associated to 

liw = [dfj,, d v \ = <9 M 7„ - <9„7 M + [7 M , 7„] , (A.6) 
will also be needed, and takes the explicit form: 

l^ij = -^(r^[Ai, Xj]) ; r M „ = -[u^u,,] - . (A.7) 
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B The heat-kernel method 



In this appendix we give a brief summary of the heat-kernel method as developed by 
Jack and Osborn j22] and provide a compilation of the results which are needed for the 
calculation. Throughout this section we work in euclidean space-time. The presentation is 
to a large extent based on [22111 and some results provided in 

Let us consider the propagator G A {x,y) associated to A, which in euclidean space-time is 
defined by: 

A x G A (x,y) = 5(x-y) . (B.l) 
In the Schwinger representation, G A (x, y) is written as an integral over the eigentime r: 

POO 

Ga(x,v) = / drp{r,e)3 A {x,y;r) . (B.2) 
Jo 

The kernel of the integral, Sa{x, V\ t ), satisfies the diffusion equation: 

d T SA(x,y,T) = -A x $ A (x, y;r) , (B.3) 

and the boundary condition S A (x, y; 0) = 5(x — y). 

SaOzj y; t) can be expanded in the Seeley-DeWitt coefficients: 

1 I I 2 °° 

Sa(x, y; r) = - — -je^^ V af(x, y)T° , (B.4) 
(4vrr)2 ^ 

and after solving the integral (jB.2J) with the regulator function p(e,r) = (47rr)^i corre- 
sponding to dimensional regularization, one gets the asymptotic formula: 

G A (x,y) = G (x-y)afi(x,y) + R 1 (x-y]Cfi)a£(x,y) 

+R 2 (x - y; cfi)a£(x, y) + R 3 (x - y; cfi)a£(x, y) 
+G A (x,y;cfi) , (B.5) 

with the coefficients: 



An 2 

Ri(x;cfi) = 2(c/x)~ £ A + iV v e 2 ' \x 

7T2 

i2, rf-i-£ 



R 2 (x;cfi) = - 2(c/i)~ £ A + A>— s— 

R 3 (x;cfi) = |x| 4 f(c/i)- £ A + iV 1 ^" 2 ," - M £ ) , (B.6) 

V 7T2 / 



where c parametrizes the renormalization prescription ( ln(c) = exp(— (ln(47r)+r'(l) + l)/2 
for MS.).). 
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Eq. ()B.5|) is only valid asymptotically for r — > and can therefore only be used to extract 
the ultraviolet behavior of the propagator. The Seeley-DeWitt coefficients are given by: 



af = —a 



ypwliw + \° 2 - \d p d p a , (B.7) 



with a and the field strength 7^ defined in Eq. (jA.5|) and (jA.6|) respectively. 
In order to calculate the divergences of the tadpole graphs like Fig. El we need the propa- 
gator with up to three derivatives, which can be calculated by use of Eq. (|B.6J) . projecting 
out the space-time independent part: 

G A (x,x) = (c/i) _£ A 2a x (x, x) + G A (x, x; c/x) , 
d*G A (x,y) ly=x = (cfiy £ A2dffl(x,y) ly=x + d*G A (x,y-,cfi) ly ^ , 

d*dlG A {x,y)\ y=x = (c^y £ k{2dldla^{x,y)\ y=x + 5^{x,x)^ 

+d x p d y l/ G A {x,y;cfi) ly=x , 

d*d x u G A (x,y) ly=x = (w)- £ A(2d*d x v a?(x,y) lv=x -6^a£(x,x) 

+dffiG A (x,y;cn) ly=x , 
dld:dy p G A (x,y) ]v=x = M^A (B.8) 



2d*d:dy p a?(x,x) 

-(5f, u d y p a£(x,x) - 5 tip d c u af{x,x) - 8 vp d%a£{x, x) ) 
+d x l d x w d y p G A (x,y;cfi) ly=x . 



For the calculation of the divergences of the functional derivative of of Zw \ Eq. (|4.5|) or 
Fig. El one needs to deal with products of propagators. After a couple of manipulations, 
the space-time dependent part of such products can be brought into the form of a sum of 
terms proportional to: 



n 



j= i * * 

Such terms can be represented by delta functions via their Fourier transforms in <i-space- 
time-dimension: 

d d z-±- = ^(-(rn-d/2)) d/2 

\z\ 2m r(m) vy ' y ' 

substituting q 2n ~ ke — > fi~ ke [—d 2 ) n 5 d (z) ; n, k G No, from which the ultraviolet divergent 
parts can easily be extracted. 
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Explicit result for 



..21 6 N 5N\ /A .. .,15 IN N\ 
+ (A X+X+X+ )(- — + J^~J + -32") + (Ax+X+MX+Xg "iv^ _ iv + T~32^ 

/a v/ v/ 1 2 3 5iV JV 2 . /A .. ..2 5 A/" . 

+(A x+ >(x + x + >( i - ^ + ^ + ^2--^-) + (Ax-Xx + x-)(-^3 - + y^) 

, A \ / 2 lx , A 2 9iV 2 7 3iV, 

+(AX_)(X + )(X_)(^4 + 2^2) + (Ax + )(x + )(x + )(- + -^4 + — - ^ - — ) 

.. .. 3 3 3A\ .. ..3 4 iV 2 . 

+(Ax_X-Xx + )(-^ - — + — ) + (Ax_X + X-)(-g + ^2 + 32") 

IN 1 2 N 2 

+(Ax + )(x_X_)(-^r - 33) + ({ A X_X_X+})(-g + + g2 ) 

+({Ax_x + })(x_)(-^ - + + ([Ax-xJXx+X^ + ^ - Yg) 

11 3 1 A" A^ 2 

+([A X _X + X+])(-g " 2jp) + (Ax + P(x + >(g - ^ + ~ Jq) 

U1W 1 N 2 9 7 3N 3N 2 

+<[A X _x + 3>(-4 + T6 } + ( A ^X + ^)(x + )(-^ + ^ - + -gy) 

/a w m/ 1 7 A^ AT 2 . /A .. uw ..1 5 1 A" 
+(Ax + Xx^X I g + ^ - j - ^ ) + (A^XxXX^X- + ^ + 8N + J ] 

/ax/ uw 1 2 ^ ^\ ,r A a iw 3 ^ A^" 2 . 

+(A« M ><x + X + « M >(-4 - + ^ + ) + ({Ax + ^X + })(-g - + -g-) 

/A X/ 2W 1 15 3iV 3iY \ /A ii\ / 1 A^" . 

+ (A X+ )( X+M 2 )(-- -_ + _ + __) + (A UllX+x+u »)(-- + Yg") 
/ a 2\ / \/ \/ 3 7 1 3AT 2w ..1 9 A" A" 2 

,.9 \/ v / 3 1 N N\ /A 2w ..3 1 JV A^ 2 

2 x, 1 3N\ /A .. 2 . . 1 19 17 N JV 2 . 
+ ( A X + - 2 X + )(-2 - -32-) + < A X + X + X- 2 X-g - g^ + - 32- + 

/a w w 2W 11 17 1 3A \ , fA 1W uw 3 JV 2 . 

+ (A X+ )( X+ )( M 2 )(- + _ + _ - -) + ({Au, x+ })( x+ un( w - — + 32) 

, r . 2 7 N 5N\ /rA u ...1 3N 3iV 2 , 

+ ( { A. 2 x + X + }>(-T6 - i6 - - fir) + ({^ X+ u» X+ })(- + ^ " ^ 

/r a 2 ,w x/ 3 1 ^ /fa a iw 9 ^ A^" 2 . 

+« A «V»<X f Xig + ^ " 7^ + 64 ) + " 32 + T } 

3 1 N N 2 IN 2 
+i([A^ x /])( x+ )(-- + _-_ + _) + i([A V /x + ])(- i + gg) 

^([A^xJXx/X-^ + + j ([ Am mX + x/])(^ + ^) 
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+ 



+ 



+ 



+ 



+ 



+ 



AV^)(x + )(^ " ^) + (Ah»„ X+ hn(-l + ^) + (A X+ ){V^>(-^ 

1 TV N 2 iV 2 
{A X+ V^»(^ " 32 + + ([AX - V^HXga ) 

9 AT N 2 N 2 
[Au mX _^ X+ ])(- + _-_) + ([A^ x+ ^ x _])(— ) 

a v/r a iw 1 3iV N /r* 2 lw 9 ^ 5Ar \ 

An M )([ X+ ^xJ)(-^ - + ([Ax_« 2 X + ])(- + - - — ) 

3 N N 2 3 1 N N 2 

lAu„ X -])(x + u»)( w - Vq - ^ ) + <[Au^]>< X+ >(- - - + - - - 

a iw 1 1 3iV ^\ /pa mw 5 5iV N 

[A^ X _^ X+ ])(-- -_-_-_) + ([Ati M X-X + «1)(-gg + ^) 

1 7 lliV 5 AT A^ 2 

[A« mX+ ])( X _^)(- + _-__) + ([An, X+ ^ X J)(-- - - + — ) 

uiw 5 -A" . ... „ ... . . 3 1 29iV AT 2 , 

ta- tt ^-ffi " 64 } + ([ A « 2 X-]><^>(32 " AN + -96" " 48 } 

[A 2 n/ 7 1 iV 3iV 2 , ... 2 ... 19 1 N 

[An X-XM-u + ^-T2 + ^2 ) + {[AX - U *+ ])( "96 + 2N 2 + 192 " 

[A 2 n/^ /pa un/29 7iV AT 2 . 

[A« 2 x + X_])( gj) + ([Ax_^X + «1)(^ + — - — ) 

iw 2w 1 1 N \ /fa 2i\/^ 3iV \ 

[Ax_X + ]>(u 2 >(- + ^ - ^) + ([Ax_X + n 2 ])(^ - -g r ) 

[Au 2 X + ]><X-X^ " ^ + f ) + KA^)(x^)(x + )(-^ + ^) 

3 1 TV iV 2 3 N N 2 

({A^x_ M })(x + )(g - ^ + Y Q ~ ig) + i({An M x/x + })(-g - ^ + -g-) 

1 7 3N 3N 2 N N 2 

(Ax + )(x + X)(g + 4/v 8 1^ } + ^ Au ^-^^^ ~ 32") 

({A^x + })(x/)(-f ) + i({Ax_Xx + })(^ " ^ + Te") 

1 N N 2 N N 2 

([A^x + n)(^ " ^) +i([A^J)(x + ^)(^ - ^) 
({Ax_,xX»(-^ + ^) + i([A W+ ]){^/)(4 + f) 

([A«^ + ^]>(-l + ^) + 1 ([Ax + ^/ i -]}(^ - g - ^) 

5 7 N N 2 

([A % x + n^1X- i g + gg) + i<[Ax + «^Xl>(i6 + 32 " Is } 
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+i<[a x+ V]> - |f ) + i([Ax + - ya - ^ 

/a v/ uw v, 1 1 3 7V , . . . . ..1 3 5N N\ 

+i(A^)(x/)(X-)(g + ^ " ^ " B ) + WVX-Xj " ^ + — - - ) 

uw 1 1 3iV N /a v/ uw 1 3 A\ 

+i<A X+/1 ><x_ U ">(- - 4^ + 16 ) + ^^-Xx+XXg + W + 16 } 
+i({A^ x _}>( x />(^ + + i({A^ x _ x /})(-i - ^) 

1 iV iV 2 5 31V ]V 2 

+i<{Ax_u„x/})(-- " is " 32) + 1 <{ A ^-X + X»(-8 " + T6 } 

+<A« # ,>( X+ V^("^ + f ) + W(X>")(-^ + ^) 

1 A/" 1 AT 2 9 AT 

+ (Au^)( x+ ^)(— + -) + (A«^ + X>(-3g " 72") + <Ax + P<u 2 >(-— - T2 ) 

+i({A^x/})(X-)(^ " ^ + ^) + (AX +M ,)(^)(-^ + ^) 

2 1 1 AT 1 AT 2 1 3 

+ (Au 2 )( X+ P(^ + — ) + <{A W+ W}>(- - — ) + <A X _>( X V>(-- - — 

93 93 /V 2 /V 2 7 N 2 

+ <{Au 2 X+ £}>(-^ + -^) + i([A X _« M x_ M ]>%) + ({An, X /X})(-^ + ^ ) 

+i([An, X / X J)(-^) + i([A^xJ)(x_ M )(^) + i{[A X _ UmX _1)(- 1 + ^) 
+i([A%_X + X])(-7 + S + (A«,)(x_^)(x_)( 



+f J/v 4 16 8 2iV 2 8iV 16 y 

1 1 N N\ .. uw 1 1 3JV, 

+(A^)( X _X_^)(-- + - - - + — ) + (A« M x_)(X-^)(- I g + ^ - 32 . 



+(A W 2 >(X-X->(-4 " i) + (AX-X_)(- 2 )(-4 + 

uw 1 1 3iV . . 2 ..5 1 N 7N\ 
+ { A x _u,)( X _un(-- + ^ - —) + (A X _u\_)(- + - + - + —) 

U TW 13 1 3N N \ /A UW 5 N ^ 

+({Au» X _u» X _})(- - - + — - -) + (Au^-X-^X^ + j) 

i r a o n/ w 1 5 iV. /fA 2 1W 1 1 N iV 2 . 
+« Au V}><X-X- I g " m ~ 33) + ({A. 2 X_X-})(3 + ^ + is - 96 J 

+{Au,)(h vp hV)(^) + <Au,)<{/^>'}>(^) + (Au,u v )(h^)(^) 

N „, , iV 2 „„,.,, iV\ 



.(A Ufl u u )(h\h^)(--) + <Au^»')(-— ) + (AurKpWWK- 



/v 72 / \ p * p - ,v 3 g/ i- -p - -- , v 144 , 

^) + (av)({^X^})(-^) + (^«V P >(~) 
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+ 



+ 



+ 



+ 
+ 
+ 
+ 
+ 
+ 
+ 



+ 



{Au, Uw h\h^})(^ - A } + {{Au^ v }){h%w){-^) + ({Au^xym^r) 

{Au,h vp })(hV){--) + ({A tt ,V}){W)(--) + ({A« 2 V^})(g + 
(A^)<[ X+ X«1>(-^)+i<[A«^ 



4g / u m -j/ w+ 48 ^ 11 y 24' 

AN 12' 



([Au^x/u^ih + i([An^ X /^])(i) + i([A V /])(u 2 }( " 



([A^X + J)(«V)(^ + ^) +i([A^ 2 ])(x/)(^f ) +i([AvVl)(^ + ^) 

1 7/V 2 1 AT 2 1 

([AuX X /})(- - +i([An 2 X+ X])(g " "ig) + i(Au,)(n^)( X /)(--) 

(An,)( X /n 2 )(- A + + i(An,)({ X _X^})(- ^ + ±A) 

({Au,u v })( X _^)( A + A) + (A«,)({ X+ A 2 }}(- A - A + A) 
(Au^Xx^X^) + i(A^x_^)(-^ - ^) + i(Au, X /)(u 2 ){- A _ 
(Au^Xu^X- A + + !(A^ X _X^)(-^ " ^) + i(Ax + ,)(nV)(^) 

( A X + X)(- 2 )(-^ - ^) + i(Ax + ^> W(- A + 
({A W ^ X /})(1 - A) + 1 ({A W ^ X _^})(-1 - A) 

(A« 2 )(x + X)(^ + ^) + ><{AW_^}>^ - ^) 

({A^ X /})< U 2 >(- A) + i({A^x_V})(-^ + ^) + i({A^x_,}>(^>(-^) 

7 A/" 2 7 83 A/" 

({A^ x _X^»(-72 + 72) + (^Xx+ti^X^ " W } 

({aa^x-H + A^) + i({a m 2 x+ x))(-^ + ^) 

<{A^ 2 x /}>(-- -—) + {A Ull ){^u v ){x + u v ){-- + -) 

Au»)(u 2 )(x + u»)(~ + + (A« M X«VXx + X-| + f ) + i({A^ 2 }><X- M >(-|A 

Am u )(y + «„//)(4 - — ) + (AuJ(y.«V)(A 7 - — ) 
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23i\r 3 iV 2 

+(A^^)( X+ ^)(- ) + (A^xX^X"^ + gg) 

„ uw 83 iV 5iV 2 . . . „ 1 137V, 

+(A W+U V>(— - - + _) + (A W „)<W>(— - —) 

+ (A^ X+ )( M V)(-^) + (Au» X+ u»)(u 2 )(^ + A) + (A V V)( x+ }(^ + ^) 

1 AT 2 1 AT 

+(A{w w y + « 2 ^})(— + — ) + (An 2 )(n 2 )(x + )(-— + — ) 
\ i J/v 144 288 ; N A /\a+/v lg lg ; 

+ <a^><^>(-- - — + ^) + <A M V)< X+ >(-- + -) 

/a 2 w 2w 1 TV . /A 2 ,. . 3 TV 1L/V 2 , 

+ AmY « 2 1 ) + Am 2 y + m 2 ) 

/v 6iV 288 ; V A+ /V 16 16 96 ; 

+(A X+ )(^^>(-f " ^) + (A X+ )(u 2 )(n 2 )(A + + (A X+ )(A 2 >(^ + ^ 
+<A X+ u M )<uV)(-^) +i(A^>(^X)( X _)(i - ^) + (A^>(^X^>(-f ) 
+i(A^)({^X x _})(^ - ^) + (A^)(nX^)(f ) + <Au M ><{V>"u'}>(^) 
+i(A v ,}(r x _}(-^) + (A^}(^>")(-^) + (A Ult u v )(h"%u"){~) 
+ (Au li u v )(h t ru")(~) + i(AnXX)(X-)(f ) +i(AV>(X-«V>(— ) 
+ ({A V ^n p ^})(^) + ({An,n^ p ^})(-1 - ^) + {{Au,u v u p h^}){-^} 
+i({A«^^})(x_)(^) + ({An^^X))(l + ^) + <{Au^>>})(i) 
+i({An^ X _^})(-i + ^) + ({A«,^})(//)(-^) 
+ ({A Ufl h/ p })(u^)(^) + i(A X _)(h fll/ uW)(^) - (AVp)(«W>(^) 

+i({A V _ u ,r})(i - + i({A^ X -^X}>(-^ + ^) 

1 iV 2 3 3A^ iV 2 

+i({A x _^^}}(- - — ) + i({A x _^ - — + — ) 

+i({A x _ V}>(«V)(-^) + 1<{A X _ V« M «"»(^ + ^) (CI) 
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(•) denotes the trace in SU(N). We calculated the generic case with N flavors in order to 
make the above counterterm also usable for the quenched case [2UI2H1- The calculation 
was performed in euclidean space-time, but above we provide the result transformed back 
to Minkowski space. Furthermore we used the notation: 

[Ax...A n ] := A x A 2 ...A n - A n A n _ x ...A 1 ; {Ai...A n } := A x A 2 ...A n + A n A n _ 1 ...A 1 . 

(2) 

The above expression corresponds only to the part of A\ , which does not involve any 
external vector sources and can contribute to processes involving four pseudo Goldstone 
fields, i.e. can be used for the (physical) K — > nit and K — > tttttt amplitudes. For ob- 
vious reasons we do not provide a minimal basis for Jlffi- However, for the part of 
shown here, all operators are linearly independent: One can use Cayley- Hamilton relations 
(CHR), the Bianchi identities and partial integrations to (potentially) reduce the number 
of operators and find a basis. Yet, the CHR are not usable since we worked in general 
SU(N), the Bianchi identities involve vector-sources which were neglected, and partial 
integrations cannot be employed since it would necessarily generate an operator involving 
A 2 

We do not to present the full result, including also the operators containing A^ and f± ' , 
here, since it is approximately a factor four times larger. It can however be obtained from 
the author. 



D Notation/Definitions 

The notation is as follows: 

U:=u :=exp(— — ) ; U — ► g R Ug l L , 

with (j), defined in Eq. (12.4)1 . representing the pseudo Goldstone bosons and the arrow 
showing the response of U to a chiral transformation (gL,9R) £ (SU(3)l, SU(3)r). 

2 In a full basis it would however be preferable to trade operators within A2 f° r operators involving 
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The building blocks used are: 



A 


:= uAqU , 




:= - ir M )M - - i/ m )m'J 


x ± 


t t i t 
:= m'xm ± MX u > 


r„ 


1 

:= 2 ( u ~~ ir /^ M + ~ X ^) M 




:= d,X+[Y^X] , 


A, 


:= V^A + ^A,^} , 




:= V^Wj, , 




:= V^u u + , 




• ^ fj,hup i 




■= V M x± - ^{Xt,M/J » 




:= V M x ± , 


X±{fi,v} 


:= {V M ,V„} X± • 


All of these transform like: 





X^h{<p,g)XtJ(<p,g) , 
where h(tp,g) is the compensating SU(N)y transformation defined: 

u -> g R uh\ip,g) = h((p,g)ug{ , 



with the exception of A^ which transforms like A^ — > g^A^g^. In addition we used notation 
that derives from the above {u 2 := u p u^, h pupa := V pV u h pa etc.). All calculations were 
performed with FORM 3.1 [2D]. 
Furthermore, we used: 

e:=A-d- iV := (4tt)- 2 ; A := — . (D.2) 
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E Expansion of the building blocks in the ^-fields 

Below we provide the expansion of the building blocks given in appendix O in terms of the 
quantum fluctuation fields £. The building blocks on the RHS are evaluated at the EOM 
(expressed by the bar). 

a = A-^[A,e]-^[[A,e],e] + ^[[[A,e],e],e] + o(e 4 ) , 
= r M +i[« M ,e]-^^e]-^pM,e],e]^] + o(e 4 ) , 

x ± = x ± -^{x T ,0-~{{x ± , 0,0 + 0,0,0 + o(£ 4 ) , 

/r = /r-^[/r^]-^[[/r,e],e] + ^[[[/r^],^,^] + o(a , 

a, = A M -i[A M ,e]-^[A M ,e],e] + ^[[[A^e]^]^] + o(e 4 ) , 

1 _ 1 
1 _ 1 1 

-^[^^^^-^^^(x-a + o^ 3 ) . (e.i) 

We use the conventions w = wexp(i£/2) , £ = A a £ a , (\ a \b) = 25 a b and the notation 
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